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The population of the Kuiper Belt within 50 AU of the Sun has likely been severely 
depleted by gravitational perturbations from the giant planets, particularly Neptune. 
The density of Kuiper Belt objects is expected to be two orders of magnitude higher 



just beyond 50 AU, where planetary perturbations are insignificant. In 1998 and 1999, 
we surveyed for Kuiper Belt Objects (KBOs) in 6 fields of the ecliptic (total sky area 



1.5 deg 2 ) to limiting magnitudes between R = 24.9 and R = 25.9. This is deep enough 
to detect KBOs of diameter > 160 km at a distance of 65 AU. We detected 24 objects. 
None of these objects, however, is beyond 53 AU. Our survey places a 95% CL upper 
limit of E < 5deg~ 2 on the surface density of KBOs larger than ~ 160 km beyond 
55 AU. This can be compared to the surface density of ~ 6deg~ 2 of > 160 km KBOs 
at distances 30-50 AU determined from this survey and previous shallower surveys. 
The mean volume density of D > 160 km KBOs in the 55-65 AU region is, at > 95% 
confidence, less than the mean density in the 30-50 AU region, and at most | of the 
mean density from 40-50 AU. Thus, a substantial density increase beyond 50 AU is 
excluded in this model-independent estimate, implying that some process or event in 
the history of the Solar System has truncated the distribution of 160-km planetesimals 
at ~ 50 AU. A dense primordial disk could be present beyond 50 AU if it contains only 
smaller objects, or is sufficiently thin and inclined to have escaped detection in our 6 
survey fields. 
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Subject headings: Kuiper Belt — minor planets — solar system: formation 

1. Introduction 

Our planetary system should not be expected to be bounded by the orbit of the ninth planet 
Pluto, as pointed out by Edgeworth (1949) and Kuiper (1951). Our current general understanding 
of the formation of our solar system is consistent with the expectation that a large population of 
small bodies - leftovers from the primordial planetesimal disk in the Solar Nebula - exists at the 
present time beyond the orbits of Neptune and Pluto where planetary accretion timescales exceed 
the current age of the solar system. Indirect evidence for such a population (now known as the 
Kuiper Belt) exists in the orbital properties of the short period comets, as has been suggested by 
several authors (Fernandez 1980; Duncan, Quinn, & Tremaine 1988). Direct observational evidence 
was first obtained with the discovery of 1992 QB1 (Jewitt and Luu 1993), and has subsequently 
grown with the detection of more than 300 objects. It is notable that all but one of the currently 
known Kuiper Belt objects (KBOs) lie within 55 AU of the Sun. 

The semimajor axis zone 30 < a < 50 AU is dynamically a very complex region. The primordial 
Kuiper Belt population here is expected to have been extensively sculpted and depleted throughout 
solar system history [see review by Malhotra, Duncan, & Levison (2000)]. Numerical simulations 
show that most primordial KBOs on orbits with a < 36 AU would have been ejected from the Kuiper 
Belt in much less than the age of the solar system (Holman & Wisdom 1993; Levison & Duncan 
1993; Duncan, Levison, & Budd 1995). In the 36 < a < 42 AU zone, the dynamical lifetime of low 
eccentricity, low inclination orbits is comparable to the age of the solar system, whereas beyond 
a ps 42 AU the dynamical lifetime is well in excess of the age of the Solar system. The orbital 
migration of the giant planets in early solar system history is expected to have rearranged the 
primordial orbital distribution in the a < 50 AU zone, by sweeping a significant fraction of KBOs 
into eccentric orbits at narrow semimajor axis zones at the locations of mean motion resonances 
with Neptune (Malhotra 1993, 1995). Furthermore, there exists a significant "scattered" KBO 
population, in very eccentric orbits with perihelia near ~ 36 AU; these KBOs are thought to have 
been formed in the 36 < a < 42 AU zone and been gravitationally scattered during a close encounter 
with Neptune (Levison & Duncan 1997). 

The current estimate of the surface density in the inner Kuiper Belt is approximately two 
orders of magnitude less than that expected from a smooth extrapolation of the surface density 
from the planetary region [cf. Weidenschilling (1977)]. It is also approximately two orders of 
magnitude less than that required for the formation of the largest observed KBOs within 10 8 yrs, 
before perturbations from the giant planets curtailed their growth (Stern 1996; Stern & Colwell 
1997b; Kenyon & Luu 1999). These arguments suggest severe depletion of the Kuiper Belt within 
50 AU. 

Beyond 50 AU, the gravitational influence of Neptune and the giant planets is insignificant. 
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A higher present-day density of objects, reflective of the primordial surface density of solids in the 
Solar Nebula, would be expected in this unperturbed distant Kuiper Belt. Previous surveys have 
detected no objects beyond 50 AU, leading several authors (Dones 1997; Chiang & Brown 1999; 
Jewitt, Luu, & Trujillo 1998) to suggest that this higher-density outer region is in fact absent, and 
that the primordial planetesimal disk must have had a cutoff near 50 AU. Gladman et al. (1998) 
calculate the fraction of detected KBOs expected to be beyond 50 AU; applying their formula to 
the number of KBOs known today also indicates that there is a deficiency of distant objects. These 
conclusions, however, depend quite sensitively upon the exponents of assumed power laws for size 
and radial distributions of the population. 

In this Letter we report the results of a survey for faint KBOs covering more than one square 
degree of the ecliptic. Most of the known KBO population lie at ~ 40 AU and have magnitudes 
R < 24. At our limiting magnitudes of R > 25, we can detect physically similar objects to a distance 
of 65 AU or greater. We compare the results of our faint survey with the detected population of 
similarly sized KBOs at smaller distances, and test whether the population density increases beyond 
50 AU. This comparison of the population of > 160 km objects is independent of any models or 
assumptions for the size distribution of the KBOs in the two regions. 

2. Search for Distant Objects 

Observations were taken using the BTC CCD Mosaic Camera (Wittman et al. 1998) at the 
Cerro Tololo Interamerican Observatory on 3 nights in May 1998 and 4 nights in May 1999. The 
BTC images 0.25 deg 2 of (non-contiguous) sky per exposure. Each field was observed for a series of 
20 to 30 eight-minute exposures spread across two consecutive nights. These images are registered 
and summed (with sigma-clipping rejection) to yield a deep template image of the fixed field. The 
template is then subtracted from each individual exposure after applying an algorithm to match 
the point spread functions, thus efficiently removing all non-moving objects from the images. These 
images are first searched for bright slow-moving objects. Then each field's images are summed with 
displacements to track a potential KBO motion, and the sum image searched for faint objects. This 
last step is repeated several thousand times to permit detection of KBOs in any bound orbit at 
distances 30-80 AU. The entire search procedure is repeated on another pair of nights six to eight 
days later. Candidate KBOs were confirmed using data taken six to eight days apart. 

A detailed description of our processing and detection methods will be published at a later 
date. 

Effective areas and limiting magnitudes are determined by inserting several thousand artificial 
KBOs into the raw data at a variety of velocities and positions. We then search for these using 
the same methods used to detect real KBOs. The total area subtended by the six search fields is 
1.5 deg 2 , but bright stars, CCD defects, etc., limit the effective area of the search to 1.3 deg 2 for 
bright KBOs. The effective search area drops at fainter magnitudes due to noise. 
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The May 1998 observations used the Kron-Cousins R filter, and we find that the effective area 
of these three fields of the search is 0.51 deg 2 at bright magnitudes, dropping by 50% at R = 24.9 to 
R = 25.4. The May 1999 observations, using a wide-band VR filter (Jewitt, Luu, & Chen 1996), 
reach 50% completeness for 25.3-25.9 R mag KBOs over effective area of 0.77 deg 2 . VR magnitudes 
are converted to R assuming a KBO color index of V — R = 0.3. Table 1 lists the fields observed 
for this survey and the total integration times, limiting magnitudes and effective areas for each. 
Several times more area is covered than for the largest previously published survey to this depth 
(Luu k Jewitt 1998). 

Twenty-three KBOs and one Centaur were discovered in the survey. The observing scheme 
produces a 10-day arc for most of these objects. Orbital parameters and their uncertainties are 
determined using the software described in Bernstein & Khushalani (2000). Although the full 
orbital parameters are poorly constrained by such a short arc, the distance is determined fairly 
accurately. This is because the object's apparent acceleration is due almost entirely to the reflex of 
the Earth's transverse acceleration, and hence the apparent acceleration is inversely proportional 
to the distance. Note that this technique requires that we not observe at opposition. Recovery 
observations have been made for 7 of the objects, and in all cases the refined orbits produce distances 
consistent with the 10-day estimates. A few of the objects were detected near the edges of the frame 
and were outside the field of view for the second (confirming) pair of nights. Distances to these 
objects with only 24-hour arcs are less certain, and dependent upon the assumption that the orbit 
is bound. Table 2 contains a list of detected objects and their heliocentric distances, and Figure 1 
plots their size vs. distance under the assumption of a 4% albedo. Distances beyond 53 AU are 
ruled out for all of these objects. Thus, over 1.3 deg 2 , we have zero detections of objects beyond 53 
AU, even though objects with diameters larger than 160 km would be visible to 65 AU over much 
of the field. 

The diameter (D [km]) of a Kuiper Belt Object can be estimated using the formula for a 
uniformly scattering sphere (Russell 1916): 

P D 2 <t) = 9 x 10 16 W^-^IO - 4 ^®"™^ (1) 
y 1AU 4 V ' 

where p is the geometric albedo, taken to be p ~ 0.04 following Jewitt, Luu, & Trujillo (1998). R 
and A are the heliocentric and geocentric distances in AU, which are approximately equal for these 
large distances; m & and tur are the apparent magnitude of the Sun and KBO, respectively; and 4> 
is a phase correction, which we ignore for these distant KBOs since the phase angle is always small. 

A quantitative characterization of the upper limit on large KBOs in the 55-65 AU range is 
as follows: At 55 AU, a magnitude of R = 25.1 is sufficient to detect D > 160 km objects, while 
at 65 AU, the necessary depth is R = 25.8. The search efficiency is dropping in this magnitude 
interval, so an exact estimate of the effective search volume would require a model for the intrinsic 
KBO distance distribution. An approximation adequate for our purposes is to take the effective 
area to be that of a D = 160 km object at the midpoint, 60 AU. This corresponds to a magnitude 
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of R = 25.5, at which the effective search area is 0.58 deg 2 . Our upper limit to the surface density 
of KBOs with D > 160 km at distances 55 < r < 65 AU is therefore 5 deg" 2 (95% CL). 

3. Comparison to the Inner Kuiper Belt 

We test the hypothesis of inner-belt depletion by comparing the volume density of KBOs inside 
50 AU to that outside 50 AU. Because we have a faint limiting magnitude and we have determined 
the distances (and hence approximate size) for all our detected objects, we can proceed by con- 
structing samples of objects of similar size in the two distance regimes, without any assumption 
about the KBO size distribution. The outer region is the 55 < R < 65 AU annulus quantified above. 

For our inner region, we select all KBOs with distances in the range 30-50 AU and diameters 
greater than 160 km. Most known KBOs within 40 AU are in mean-motion resonance orbits or 
have likely been perturbed by a Neptune encounter, affecting the surface density in an uncertain 
manner. For this reason, we will also compare our results to an inner region of 40-50 AU. 

Our survey yields 8 KBOs which have a distance between 30-50 AU and diameter greater than 
160 km. These are all brighter than R = 24.1, although at 50 AU a 160 km KBO is only R = 24.7. 

Denoting the mean volume density of D > 160 km KBOs in the inner and outer regions as 
n\ and n2, respectively, we are interested in an upper bound to the density ratio / = U2jn\. We 
assume that A^i objects were detected in the inner region over an effective volume Vi, and N2 = 
objects were detected in the effective volume V2 of the outer region. A Bayesian analysis with 
uniform prior on n\ yields the following simple formula for the probability of / being above some 
value: 

p(>f) = (i + m/v 1 )- N > (2) 

The effective volume for the inner region is simply that subtended by our survey's 1.3 deg 2 effective 
area over the 30 < R < 50 AU distance range, since the inner-region KBOs are well above our 
detection threshold. The effective volume of the outer region is given by the 0.58 deg 2 effective 
area through the 55-65 AU depth. The ratio of volumes between the 55-65 AU sample and the 
30-50 AU sample, V2/V1, is then 0.49. Equation (2) then gives / < 0.92 at 95% confidence. If only 
the 7 KBOs in the inner region from 40-50 AU are considered, then we find / < 0.67 (95% CL). 
We can thus assert in a model-independent way that the density of KBOs in the 55-65 AU range 
is at least | lower than in the 40-50 AU range. 

We could estimate the inner-region density from other surveys with brighter magnitude limits 
than ours, though the comparison with our own brighter objects has the advantage of cancelling 
any dependence on ecliptic latitude or longitude which may be present in the KBO density. There 
are very few published wide-area KBO surveys to serve as comparison samples. Jewitt, Luu, & 
Chen (1996, JLC) and Jewitt, Luu, & Trujillo (1998, JLT) are the two largest. The JLC survey 
covers 8.3 deg 2 to a limiting magnitude of R = 24.2. At this depth, KBOs must have diameter 
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D > 200 km, to be visible all the way to 50 AU. Out of 15 KBOs found by JLC, 12 have distances 
within 30-50 AU and D > 200 km. Using JLC's effective areas quoted at R = 23.2 and R = 24.2, 
a rough estimate for the effective area for all 12 KBOs is 6.1 deg 2 . Our survey will detect 200 km 
objects to a distance of 67 AU, if the same R = 25.5 limit is used, so we may conservatively take 
the outer region in this case to be 55-67 AU across the 0.58 deg 2 effective area for R < 25.5. This 
then yields a Bayesian limit of / < 2.17 (95% CL) for the ratio of outer to inner volume densities 
of D > 200 km objects. 

The JLT survey, with a much brighter limiting magnitude of R = 22.5 over 51.5 deg 2 , detects 
13 KBOs. Of these, 12 are between 30-50 AU, but only 6 are larger than D = 320 km, the smallest 
size that can be seen at 50 AU in this survey. These 6 KBOs are all inside 45 AU, although we 
will still consider the inner region 30-50 AU. Objects with D > 320 km are detectable to 85 AU at 
R < 25.5, so we take the outer region to be 55-85 AU. Using the same procedures as for the JLC 
survey, we calculate / < 4.9. 

We should at some point expect to find objects beyond 55 AU since the orbits of the known 
scattered-disk members will carry them well beyond this point. The sole known KBO that is likely 
beyond 55 AU (1999 DGs) could be a scattered-disk member; its 1-night arc is insufficient to 
determine an orbit. Seven objects with thr > 25 were discovered in very deep Keck pencil-beam 
surveys (Luu & Jewitt 1998; Chiang & Brown 1999); all have motions over arcs of a few hours 
that are consistent with distances < 50 AU. 

Our failure to detect any scattered-disk objects beyond 50 AU does not invalidate our argument, 
since we claim only a fair sample of the population, not an exhaustive survey. We can test the 
fair-sample hypothesis by examining our detected population of D > 160 km objects; we claim 
that this population is complete to about 65 AU, and that majority of the population is within 
55 AU, i.e. that we are capable of detecting most of the objects anywhere in their orbits. There 
should consequently be no bias toward finding objects at perihelion. Examining the 9 objects with 
D > 160 km in Table 2, we see that 4 are closer to aphelion than perihelion (1998 KYq±, KG62, 
KR.65, and 1999 KRig), 1 is nearer perihelion (1998 KS65), and 4 have uncertain a (1999 JA132, 
JB132, JD132, and JF132). Our objects' distance distribution is thus qualitatively consistent with 
an unbiased sampling of their own orbits. A quantitative comparison (KS test) also demonstrates 
internal consistency. While this is a weak test, it is one which the full population of known KBOs 
fails miserably — there is an extremely strong bias toward objects near perihelion. Indeed our 1999 
KRig is the only known KBO of significant orbital eccentricity (e > 0.1) to be discovered near 
aphelion. 

Our lack of detections beyond 55 AU is consistent with a density beyond 55 AU that is at 
most similar to that within 50 AU, and it is likely we have observed an under-density in the outer 
Kuiper Belt. The measurements are not consistent with a large increase in the surface density 
beyond 50 AU. 
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4. Conclusions 

If planetary perturbations are solely responsible for the structure of the Kuiper Belt, a dense 
primordial disk would be expected beyond ~ 50 AU where these perturbations are insignificant. 
Our survey could have detected such a disk but did not. 

There are several possible explanations for this non-detection. 

1. KBOs in the outer Kuiper Belt could be fainter than expected, owing to a lower albedo or 
much redder color or much smaller sizes. To explain our survey results, this would require a 
change in the physical properties of KBOs beyond 50 AU. 

2. The outer Kuiper Belt could have been dynamically excited early in the history of the solar 
system by a stellar encounter (Ida et al. 2000) most likely during the Sun's residence in 
its birth cluster (Adams & Laughlin 2000); by perturbations from large Neptune-scattered 
planetesimals (Petit et al. 1999) or by proto-planetary cores (Thommes et al. 1999). Such 
excitation could have increased the orbital inclinations and eccentricities of the objects, low- 
ering the apparent surface density on the sky and possibly decreasing the maximum size of 
objects through the cessation or reversal of the accretion process. More extreme excitation 
could have stripped away the outer Kuiper Belt. 

3. The disk is actually present, but dynamically cold (Hahn 2000). It could have escaped our 
survey by simply not intersecting our fields. This is entirely possible, if this disk is inclined 
to the ecliptic plane by as little as 1°. The invariable plane (the angular momentum plane of 
the solar system), inclined to the ecliptic by 1.5°, would be a likely candidate for the location 
of such a cold, dense disk. Our Field G is only 0?6 from the invariable plane, so the scale 
height of the cold disk would have to be <C 1° to have escaped detection. This implies a very 
thin, extremely dense disk. Absence of > 50 AU detections in future deep survey fields would 
eliminate the possibility of a cold disk. 
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grant #NAG5-6886. 



REFERENCES 

Adams, F., & Laughlin, G. 2000, Icarus (submitted). 



-8- 



Bernstein, G. M., k Khushalani, B. 2000, AJ (in press) 
Chiang, E. I., & Brown, M. E. 1999, ApJ, 118, 1422 

Dones, L. 1997, in ASP Conf. Proc. 122, From Stardust to Planetesimals, ed. Y. J. Pendleton & 
A. G. G. M. Tielens (San Francisco: ASP), 347 

Duncan, M., Quinn, T., & Tremaine, S. 1988, ApJ, 328, L69 

Duncan, M.J., Levison, & H.F., Budd, S.T. 1995, AJ, 110, 3073 

Edgeworth, K. E. 1949, MNRAS, 109, 600 

Fernandez, J. A. 1980, MNRAS, 192, 481 

Gladman, B., Kavelaars, J. J., Nicholson, P. D., Loredo, T. J., & Burns, J. A. (1998) AJ, 116, 2042 

Hahn, J.M. 2000, Lunar and Planetary Science XXXI 

Holman, M.J. & Wisdom, J. 1993, AJ, 105, 1987 

Ida, S., Larwood, J., & Burkert, A. 2000, ApJ, 528, 351 

Jewitt, D., and Luu, J. 1993, Nature, 362, 730 

Jewitt, D., Luu, J.X., & Chen, J. 1996, AJ, 112, 1225 

Jewitt, D., Luu, J.X., & Trujillo 1998, AJ, 115, 2125 

Kenyon, S.J. & Luu, J.X. 1999, AJ, 118, 1101 

Kuiper, G. P. 1951, in Astrophysics: A Topical Symposium, edited by J. A. Hynek (McGraw-Hill, 
New York) 357 

Levison, H.F. & Duncan, M.J. 1993, ApJ, 406, L35 

Levison, H.F. & Duncan, M.J. 1997, Science, 277, 197 

Luu, J.X. & Jewitt, D. 1998, ApJ, 502, L91 

Malhotra, R. 1993, Nature, 365, 819 

Malhotra, R. 1995, AJ, 110, 420 

Malhotra, R., Duncan, M., & Levison, H. 2000, in Protostars and Planets IV (Tucson: University 
of Arizona Press), 1231 

Petit, J., Morbidelli, A., & Valsecchi, G.B. 1999, Icarus, 141, 367 

Stern, S.A. 1996, AJ, 112, 1203 



-9- 



Stern, S.A. & Colwell, J.E. 1997a, ApJ, 490, 879 
Stern, S.A. & Colwell, J.E. 1997b, AJ, 114, 841 

Thommes, E.W., Duncan, M.J., & Levison, H.F. 1999, Nature, 402, 635 
Russell, H.N. 1916, ApJ, 43, 173 
Weidenschilling, S.J. 1977, Ap&SS, 51, 153 



This preprint was prepared with the AAS IATjrjX macros v5.0. 



-10- 




Fig. 1. — Diameter vs distance for the KBOs discovered in this survey. Squares are objects dis- 
covered in 1998, triangles are from 1999. The two curves denote lines of constant magnitude, 
corresponding to the 50% completeness limits of the overall and faintest survey fields. The two 
rectangular boxes contain denote the "inner" and "outer" regions of the Kuiper Belt described 
in the text. The absence of any objects in the outer sample region limits the volume density of 
> 160 km diameter objects at 55-65 AU to be below that at 40-50 AU at 95% CL. 



Table 1. Field Information 
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Observations marked with asterisks were searched for KBOs; unmarked observations were used only for (p)recovery. 
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Table 2. Objects Discovered 
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F 


9 d 


24.1 


32.4 ±5.0 


0.20 ±0.22 


29.8 ±6.4 


39.1 ± 1.4 


129 


1999 JF132 


F 


9 d 


24.0 






1.6 ±0.4 


43.1 ±2.1 


164 


1999 JH132 


F 


9 d 


25.5 






0.6 ±0.2 


41.1 ±2.4 


75 


1999 JJ132 


F 


9 d 


25.5 






3.2 ± 1.2 


50.1 ±2.8 


111 


1999 JK132 


F 


9 d 


24.8 






16.0 ±5.9 


39.0 ±2.5 


93 


1999 KT16 


F 


1 d 


25.1 






8.5 ±3.7 


46.3 ±3.1 


114 


1999 KK17 


F 


9 d 


25.0 






8.7 ±16. 


49.8 ± 7.5 


139 


1999 KL17 


G 


90 d 


25.2 


46.2 ±0.2 


0.03 ±0.17 


2.8 ±0.0 


47.6 ±0.0 


115 


1999 KR18 


G 


89 d 


24.9 


43.3 ±0.3 


0.21 ±0.03 


0.6 ±0.0 


52.6 ±0.1 


162 


G3-01 b 


G 


9 d 


25.6 


39.9 ±2.7 


0.16 ±0.05 


1.6 ±0.1 


33.4 ±1.4 


47 



a No data are given when the arc is too short to provide meaningful constraint. 

b Objects E2-01 and G3-01 were not reported to the MPC due to insufficient S/N on the recovery obser- 
vations. 

c Diameters assume albedo of 0.04. 



